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Abstract 

It is well known that an irreducible algebraic curve is rational (i.e. parametric) 
if and only if its genus is zero. In this paper, given a tolerance e > and an 
e-irreducible algebraic affine plane curve C of proper degree d, we introduce the 
notion of e-rationality, and we provide an algorithm to parametrize approximately 
affine e-rational plane curves, without exact singularities at infinity, by means 
of linear systems of (d — 2)-degree curves. The algorithm outputs a rational 
parametrization of a rational curve C of degree at most d which has the same 
points at infinity as C. Moreover, although we do not provide a theoretical 
analysis, our empirical analysis shows that C and C are close in practice. 
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Introduction 



Let 0* be an algebraic or geometric object that satisfies a property ^3 that implies 
the existence of certain associated objects O*; for instance, 0* might be a polynomial, 
*}3 the fact of being reducible and 0* the irreducible factors. Computer algebra tech- 
niques provide, for a wide class of situations, algorithms to check *}3, and to compute 
exactly the associated objects 0*. However, in many practical applications, we receive 
a perturbation of 0*, where does not hold anymore neither the associated objects 
0* exist. The problem, then, consists in computing a new object 0, close to and 
satisfying <}3, as well as the associated objects 0i to 0. We call approximate to an 
algorithm solving a problem of the above type. Here, the notion of "closeness" depends 
in general on the particular problem that one is solving. 

One can find in the literature approximate algorithms for computing gcds (see [3], 
[Z], [IB]), factoring polynomials (see [I], [10], [15], [21]), etc. For algebraic varieties there 
also exist approximate solutions: see [5], [6] for the implicitization problem, in [Hj the 
numerical condition of implicitly given algebraic curves and surfaces has been analyzed, 
and see [2], [9], [12], [17], [18] where the parametrization questions are treated. 

In this paper we consider the approximate parametrization problem for affine plane 
algebraic curves. That is, with the above terminology, 0* is an affine plane curve, ^3 is 
the fact of being rational, and 0* is a rational parametrization of 0*. So, the problem 
is stated as follows: we are given an affine curve (say that it is a perturbation of a 
rational curve) and we want to compute a rational parametrization of a rational affine 
curve near it; where we use the notion of "vecinity" introduced in |17j . 

In [17] and [18] the approximate parametrization problem is solved for the special 
case of affine plane curves and affine surfaces being a perturbation of a monomial 
curve and surface, respectively. In both papers, the basic tool is the use of e-points 
(see also [19]). More precisely, given a tolerance e > 0, in [T7J, the parametrization 
problem is solved for the case of affine plane curves having an e-singularity of maximum 
multiplicity, and in |18J the problem is solved for affine surfaces having also an e- 
singularity of maximum multiplicity. The basic idea was to use a pencil of lines through 
the e-singularity and, hence, it was solved working as in the exact case for monomial 
varieties. 

In this paper, we generalize the ideas in [T7J to the case of affine plane curves with- 
out singularities at infinity. For this purpose, the first obstacle is to associate suitably 
the different e-singularities. This leads to the notion of cluster. Then, we introduce 
the notion of (affine) e-rationality, and we provide an algorithm to parametrize ap- 
proximately e-rational curves without exact singularities at infinity. The idea of the 
algorithm is to work with linear systems of curves of degree d — 2, where d is the degree 
of the input curve. This system plays the role of the linear system of adjoint curves 
in the exact parametrization algorithm. In addition, we prove that the degree of the 



2 



output rational curve is bounded by the degree of the input one, and that both curves 
have the same points at infinity. Differently to [17] we do not provide a theoretical 
analysis of the error (i.e. on the closeness of input and output). However, our empirical 
analysis shows that the curves are in practice near, and it allows us to think about a 
theoretical treatment of this fact future project. 

The paper is structured as follows. In Section 1 we recall the main notions and 
properties on e-singularities. Section 2 is devoted to recall the main ideas of the exact 
parametrization algorithm for curves. In Section 3 we develop the idea of cluster 
and we introduce the notion of e-rationality. In Section 4 we derive the approximate 
algorithm, as well as the main properties of the output curve. In Section 5 we illustrate 
the algorithm by some example, and in Section 6 we analyze empirically the error. 

Throughout this paper, we use the following terminology. || • || and 1 1 - 1 1 22 denote 
the polynomial oo-norm and the usual unitary norm in C 2 , respectively. | • | denotes 
the module in C. The partial derivatives of a polynomial g G C[x, y] are denoted 
by 9 V := Qi x Qj y where = G N 2 ; we assume that g° = g. Moreover, for 

~v = G N 2 , \~v\ =i + j. Also, et = (1,0), and = (0, 1). 

In addition, we use the following general assumptions. A tolerance e is fixed 
such that < e < 1. C is an affine real plane algebraic curve over C of proper degree 
d > (see Def. 11.11) . without (exact) singularities at infinity, not passing through 
(1 : : 0), (0 : 1 : 0), and defined by an e-irreducible polynomial f(x,y) G K[x, y] ; that 
is / can not be expressed as f(x,y) = g(x,y)h(x,y) +8(x,y) where h,g,£ G C[x,y] 
and ||£(a;, ?/)|| < e||/(x, y)\\ (see [3], [I3])- We denote by C h the projective closure of C. 

Let us mention that, although we require that C is real, the results in this paper are 
also valid for non-real plane algebraic curves. In addition, the condition (1 : : 0), (0 : 
1 : 0) G" C h can be avoided by performing a suitable affine orthogonal linear change of 
coordinates. The requirement on the smoothness of C h at infinity, might be avoided by 
performing a suitable projective linear change of coordinates. However, differently to 
affine orthogonal linear changes, in general, projective changes of coordinates do not 
preserve properly the closeness between the input and output curves. 

1 Preliminaries on e-points 

Our fundamental technique to deal with the approximate parametrization problem is 
the use of e-points. The notion of e-point of an algebraic variety was introduced by 
the authors (see [IT], [IB], [19]) as a generalization of the notion of approximate root of 
a univariate polynomial. In this section, we briefly summarize some previous notions 
introduced in [T7] and [IB], and geometric properties obtained in [Tj5]. We start with 
the notion of proper degree. 

Definition 1.1. We say that a polynomial g G C[x, y] has proper degree £ if the total 
degree of g is £, and 3 v G N 2 , with | v | = £, such that \g v \ > e\\g\\. 
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We say that an algebraic plane curve has proper degree l if its defining polynomial 
has proper degree £. □ 

The notion of e-point is as follows. 

Definition 1.2. P G C 2 is an e-(affine) point of C if |/(P)| < e||/||. □ 

In this situation, we introduce the notion of e-singularity, pure e-singularity, and 
e-ramification point. 

Definition 1.3. Let P G C 2 be an e-point of C. 

(i) The e-multiplicity of P on C (we denote it by mult e (P, C)) is the smallest natural 
number r G N satisfying that 

(1) G N 2 , such that < |lf | < r - 1, it holds that |/*(P)| < e||/||, 

(2) 3 V G N 2 , with |"y | = r, such that \f~*{P)\ > e\\f\\. 

(ii) P is an e-(affine) simple point of C if mult e (P, C) = 1; otherwise, P is an e-(affine) 
singularity of C. 

(iii) P is a /c-pure e-singularity of C, with fc G {1,2}, if mult e (P, C) > 1 and 
|r ulte(P ' C) ^(P)|>e||/||. 

(iv) P is an e-(affine) ramification point of C if mult e (P, C) = 1, and either |/ ei (P)| < 
e||/||or|/-(P)|<e||/||. □ 

Note that, since C has proper degree, < mult(P, C) < mult e (P, C) < deg(C), 
where mult (P,C) denotes the "exact" multiplicity of P on C. For instance, the origin 
has exact multiplicity 1, and e-multiplicity 2, on the curve defined by |x + x 3 + y 2 . 
In the exact case, if C is irreducible, mult (P,C) < deg(C). Thus one may expect that 
in the approximate case, if C is e- irreducible, then mult e (P, C) < deg(C). Although 
this is the case in all the examples we have tried, we have not been able to prove 
it. So in this paper, when computing e-multiplicities, we also consider the possibility 
mult e (P,C) = deg(C). 

The following lemma is a direct generalization of Lemma 3 in [17] . 

Lemma 1.4. Let P G C 2 be an e-point ofC. There exists 5 > such that every Q G C 2 , 
satisfying that \\P — QW2 < 5, is an e-point of C with mu\t e (Q,C) > mult e (P, C). 
Proof. Simply observe that the reasoning of Lemma 3 in [TTj is also valid over C. □ 

The following example shows that, in Lemma 11.41 the e-multiplicity of Q can be 
strictly bigger than mult e (P, C). 



4 



Example 1.5. Let C be defined by f(x, y) = x 3 y + y 3 x + x 3 + |x 2 + ey + |; note that 
ll/H = 1. For P = (0,0), one has 

/(P) = |, /^(P) = 0, r(P) = e. 

So, mult e (P, C) = 1. Now, we consider the sequence of points {P n = (— ^,0)} n >i. 
Then, 

2 2n z n 6 n z n n 6 

n rr 

So, for n sufficiently large, mult e (P n ,C) = 3. □ 

Definition 1.6. Let P be an e-point of C and r = mult e (P, C). If P is fc-pure, with 
6 {1, 2}, we define the £;-weight of P as 



weight fc (P) = max 

t=0,...,r — 1 



r! ■ /*'^(P) 



i 

r — i 



We define the weight of P, denoted by weight(P), as max{weight 1 (P), weight 2 ((P)}, if 
P is pure in both directions, and as the corresponding fc-weight otherwise. □ 

The following two rational functions were introduced in [22], and will play an im- 
portant role in this development: 

16x \ m . . 1 x(l - 9x) 32x 2 
TZ- m (x)=2x[ z — + - — ], Tl out (x) 



l + 3x (l + 3x)V' 2 2(1 + 3x) (l + 3x) 3 ' 

Furthermore, these two rational functions give information on how close an e-point is 
to an exact point of the curve C (see Theorem 2 and Corollary 3 in |19j). 



2 Preliminaries on Symbolic Parametrization 

In this section, we briefly recall the symbolic parametrization algorithm for rational 
plane algebraic curves of degree d > 2 (note that lines and conies can be trivially 
parametrized by lines) based on (d — 2) adjoint curves; for further details see [23]. For 
this purpose, throughout this section we assume that C is rational (i.e. its genus is 
zero). In addition, taking into account our requirements in Section 131 and for simplicity 
sake, we assume in this section that all singularities are affine and ordinary. Again, for 
a complete description see [24] . 

The idea is to use a linear system of curves such that for almost every curve in this 
system, all its intersections with C h , except one, are predetermined; recall that C h is the 
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projective closure of C. Moreover, the set of all these intersection points is the same one 
for every curve in the system, and the points in this set are called the "base points" . 
Thus, if one computes the intersection points of C h with a generic representative of the 
system, the expression of the unknown intersection point gives the parametrization of 
the curve in terms of the parameter defining the linear system. 

More precisely, let Hd-2 be the linear system of adjoint curves to C h of degree d — 2. 
That is, Tid-2 is the linear system of curves of degree d — 2 having each r-fold of C h as a 
base point of multiplicity r — 1; i.e. as a point of multiplicity at least r — 1. In particular 
it implies that the multiplicity of intersection of a curve in Tid-2 and C h at a base point 
of multiplicity r — 1 is at least r(r — 1). Thus, using that the genus of C is zero, and 
taking into account Bezout's Theorem, one deduces that d — 2 intersections of C h and a 
generic element in 7id-2 are not predetermined. In this situation, one may take (d — 3) 
simple points on C h , and determine the 1-dimensional linear subsystem 7~C d _ 2 of H.d-2 
obtained when these simple points are required to be base points of multiplicity 1. 
In this way, the number of predetermined intersections (counted with multiplicity) is 
(d — l)(d — 2) + (d — 3), i.e. only one intersection point is missing. Thus, computing 
this free intersection one finds a rational parametrization of C h . Summarizing these 
ideas one has the following process: 

(1) Compute the singularities of C h as well as their multiplicities (recall that we have 
assumed that all singularities are affine and ordinary). 

(2) Determine the linear system H.d-2 of adjoint curves of degree (d — 2) to C h . 

(3) Compute d — 3 different simple points on C h . 

(4) Determine the linear subsystem H d _ 2 of TLd-2 by requiring that every simple 
points in Step 3 is a base point of multiplicity one. 

(5) Compute the free intersection point of / H d _ 2 and C h . 

Let us make a comment on how to computationally perform the steps in the above 
process. Step (1) can be performed, for instance, using resultants. In Step (2), one 
considers a homogeneous polynomial H(x,y,z) of degree (d — 2) with undetermined 
coefficients. Now, for each singular point P of multiplicity r one requires that H and 
all its partial derivatives till order (r — 1) vanish at P. This generates a linear system of 
equations in the undetermined coefficients of H. Solving it, and substituting in H, we 
get the defining polynomial of Tid-2-, let us call it again H. Step (3) may be performed 
by intersecting C h with lines (see [24J for advanced approaches); although it is not 
necessary, looking for the parallelism with the reasoning in Section HI we take affine 
simple points. Step (4) can be approached as Step (2), i.e. requiring that H vanishes 
at each simple point, solving the provided linear system and substituting the solution 
in H; let H*(t, x, y, z) be the defining polynomial of H d _ 2 (note that dim(7i^_ 2 ) = 1) . 
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Finally, let us deal with Step (5). For this purpose, let {Qi := (q^i : : l)}j=i,..., s be 
the singularities and r, the multiplicity of Qi. Also, let {Pi := (p^i : p^2 : l)}i,...,d-3 be 
the simple points determined in Step (3). Then, the free intersection point is obtained 
by computing the primitive part, w.r.t. t, of the resultants of H*(t,x,y, 1) and f(x,y) 
with respect to x and y, respectively. Indeed, it holds that (see [20] ) 

S 1 (x, t) = Res y {H*(t, x, y, 1), f(x, y)) = Yl° =1 (x - q iA y*^ ntf(x - Pi,i)M x (x, t), 
S 2 (y,t) = Res x (H*(t,x,y, 1), f(x,y)) = UUiv ~ QiflY*^ Y^Zliv -Pw) M 2 (y,t), 

where deg x (Mi) = deg y (M 2 ) = 1. Therefore, the parametrization is the solution in 
{x,y} of{M 1 (x,t) = 0,M 2 (y,t) = 0}. 

3 e- Rational Curves 

In this section we introduce the notion of e-rationality of a plane algebraic curve. This 
notion plays the corresponding role in the approximate frame that the rationality does 
for exact algebraic curves. We will have two main difficulties. On one hand, computing 
the e-multiplicity and on the other, differently to the exact case, we will have in general 
more e-singularities than expected, and we will need to associate them; we will solve 
this last difficulty introducing a suitable concept of cluster. 

We first need to determine the e-singularities. To check the existence and perform 
the actual computation of the e-singularities, one has to solve the system of algebraic 
equations 

{f t (x,y) = 0, f*(x,y) = 0, f*{x, y) = 0}, 

under fixed precision e||/||. This can be done, for instance, by applying root finding 
techniques (see [3], [11], [13]). Note that since / is e- irreducible then it is irreducible, 
and hence the above system has finitely many solutions. Let S\ be the set of solutions. 
One may accelerate the computation by working (if possible) with two co-prime poly- 
nomials, instead of three, to get a finite super-set of the set of solutions, from where 
the e-singularities are detected afterwards. 

Now, for P G Si, we want to compute mult e (P, C). This can be obviously done 
by substituting P at the corresponding partial derivatives and checking the conditions 
in Def. 11.31 (1). Seemingly, there is no difficulty on that. However, in Lemma 11.41 
we have seen that for each e-point P of e-multiplicity r there exists an open disk U 
centered at P such that if Q G U, then Q is an e-point of e-multiplicity at least r. 
So, an small perturbation of P may produce an incorrect answer for the e-multiplicity; 
see, for instance, Example 11.51 We are indeed interested in assigning the maximum 
possible e-multiplicity to the e-point. The proof of Lemma 3 in [17], and hence of 
Lemma 11.41 shows how to detect the radius of one of these open disks, so one may try 
to estimate the maximum e-multiplicity at the disk. Nevertheless, in practice, this is 
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unfeasible. Instead, we propose a different strategy that, although it does not ensure 
the achievement of the maximum, in practical examples turns to work efficiently. 

More precisely, for each k G {2, . . . , d — 1}, we take ui, . . . , u s G N 2 , with 2 < s < 
k + 1 (in practice s = 2) such that for all i, \ui\ = k and gcd(/ M1 , . . . , f Us ) = 1, and we 
solve {f ui = 0, . . . , f Us = 0}, under fixed precision e||/||. Let A k be the set of solutions. 
Then, for k G {2, . . . , d — 1} we consider the set (note that Si is defined above) 

S k = {P G A k /\f{P)\ < e||/|| \/w G N 2 with |^| < k}. 

If for a given k and for all s it holds that gcd(/ Ml , . . . , f Us ) ^ 1, we take S k = 0. Finally 
we consider the set 

5 = |J S k . 

k=l 

It is clear that in general we introduce additional e-singularities, and we will have to 
generate a process (the cluster construction) to identify them. Nevertheless, each new 
e-singularity, after identification, will increase the e- multiplicity of the original one. 

Definition 3.1. The set S, introduced above, is called the e-(affine)-singular locus of 

C. We denote it by Sing e (C). □ 

Example 3.2. Let us take e = 0.001 in Example II .51 The e-singular locus of this curve 
is Sing e (C) = Si U S 2 U S 3 where 



<Si = {Pi 
Pi 
Pz 
Pi 
P 5 
P 6 
Pi 

5 2 = {Ps 

5 3 = 0. 



0.02131893405 + 0.009609927603z, 0.02442855631 + 0.1171004584z), 
0.004713033954 + 0.02355323617i, -0.07491796596 - 0.09032199938i) 
-0.01424770212 + 0.01818884517i, 0.1084633939 + 0.05315246871?), 
-0.02443272919, -0.1159479025), 

-0.01424770212 - 0.01818884517i, 0.1084633939 - 0.05315246871i), 
0.004713033954 - 0.02355323617i, -0.07491796596 + 0.09032199938?) 
0.02131893405 - 0.009609927603i, 0.02442855631 - 0.1171004584i)}, 
-0.0001666666667,0)}. 



Moreover, mult e (Pi) = ■ • • = mult e (P7) = 2 but mult^Pg) = 3. Note that considering 
only Si we would have not found a point with e-multiplicity 3. □ 

As we could check in the previous example, the difficulty appears when observing 
that we may have two (in general more than two) e-singularities P and Q that are 
very " close" , and somehow we need to identify them. To approach this, we introduce 
the notion of cluster of e-singularities. Intuitively, two e-singularities P and Q of C are 
in the same cluster, if the disks centered at P and Q (of certain radius) are a small 
vibration of each other. The radius and the vibration are measured by means of the 
value of the function TZ OVLt at the weight and the tolerance, respectively (see Section 
[TJ). Since the notion of weight requires that the e-singularities are pure, for non-pure 
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e-singularities we will take radius zero. More precisely, we introduce the following 
definition. 

Definition 3.3. Let P be an e-point of C. We define its radius, and we denote it by 
radius(P), as TZ out (weight (P)) if P is pure and zero otherwise. □ 

Definition 3.4. Let A be a finite set of e-points of C. For P G A we define the cluster 
of P w.r.t. A as the set of all points Q e A such that at least one of the following 
conditions is verified: 

(1) ||P - Q|| 2 + | radius(P) - radius(Q) |< K out {e), 

(2) there exists P' G A such that ||P' — P||2+|radius(P') — radius(P)| < 7£ ou t(e) and 
||P' - Q|| 2 +|radius(P') - radius(Q)| < K out (e). 

We say that R is a candidate to be the representative of a cluster, if R is a point of 
the cluster of maximum e-multiplicity. We say that R is a representative of a cluster 

if it is a candidate and |/(P)| < \f(Q)\ for all the other candidates Q. We define the 
e-multiplicity of the cluster as the e-multiplicity of any of its representatives. 

We denote a cluster by Cluster r (R, A), where r is the e-multiplicity and R a rep- 
resentative, and by Cluster r (R) when A = Sing e (C). □ 

Now, we are ready to introduce the notion of e-rationality. 

Definition 3.5. If {Cluster ri (Pi)}i = i^.^ s is the cluster decomposition of Sing e (C), we 
say that C is e-(affine) rational if (d — l)(d — 2) — J2i=i r i( r i — 1) = 0. □ 

Remark 3.6. Note that in the previous theoretical development we have not consid- 
ered singularities (neither e-singularities) at infinity. We leave this extension of the 
concept of e-rationality for further research. □ 

If we apply the previous ideas to Example 13.21 (see also Example 11.51) . with e = 
0.001 we get that the 8 points of Sing e (C) belong to the same cluster. So, the cluster 
decomposition is {Cluster 3 (P 8 ) = {Pi, . . . , P 8 }}. Therefore, C is e-rational; indeed, 
it is e-monomial, and thus parametrizable with the techniques in [17]. We finish the 
section with a more general example. 

Example 3.7. Let us consider e = 0.005 and the curve C of proper degree 5 defined 

by the polynomial (see Fig{T]): 

f(x,y) = -2.199771784x 2 - 0.2197717843x 4 y - 0.9016804979:rV + 1.858817427x 3 - 
1.8916804982/ 4 +0.9899999999xy 3 +0.9899999999x 2 y-1.055726141x 2 y 3 +0.3409543568y 2 + 
0.9899999999x 4 + 0.9899999999xy 4 + 0.9899999999y 3 - 0. 1869087137x 5 + 5.235497925xy 2 - 
1.770497925xV + 1.45213693x 3 y - 0.1440456432xy - 0.52786307y 5 + 0.01. 
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The e-singular locus is Sing e (C) = Si U <S 2 U S3, where 



Si 



s 2 
s 3 



{P 1 = (-0.9956027274 + 0.0004067223817*, 0.001447687187+ 0.9982777543z), 
P 2 = (1.011706789 - 0.1320874194z, -1.008532436 + 0.06832949372z), 
p 3 = (1.007458642, -1.044045331), P 4 = (0.9909273695, -0.9540334161), 
P 5 = (1.011706789 + 0.1320874194z, -1.008532436 - 0.06832949372*), 
P 6 = (-0.9956027274 - 0.0004067223817z, 0.001447687187 - 0.9982777543i), 
A= (0,0)}, 

{P 8 = (1.000000001, -1.)}, 



Moreover, mult,, (Pi) = 
mult e (P 5 ) = mult e (P 6 ) 
(see Fig. [1]): 



mult e (P 2 ) = mult e (P 7 ) = 2, and mult e (P 3 ) = mult e (P 4 ) = 
= mult e (P 8 ) = 3. Furthermore, the cluster decomposition is 



Muster 2 (Pi) 
Cluster 2 (P 2 ) 
Cluster 2 (P 7 ) 
Cluster 3 (P 8 ) 



{Pi}, 
{P2}, 
{P7}, 

{P3, P^ P5, P&, P&}- 



Thus, C is e-rational. 





Figure 1: Left: Clusters. Right: Curve C 



□ 



4 Approximate Parametrization Algorithm 

In this section, we present our approximate parametrization algorithm. For this pur- 
pose, we assume that C is e-rational of proper degree d > 2 (note that for d — 1 the 
problem is trivial, and for d = 2 one can apply the algorithm in [17]), and that 



{Cluster n (Qi)} 



i=l,...,s. 



where Qi :- 



%2 ■ 1) 



is the cluster decomposition of Sing e (C). Furthermore, if possible, i.e. when there exists 
a real representative of the cluster, we take Qi real. 
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In this situation, we adapt the algorithm in Section [2] as follows. Let C h be the 

projective closure of C. We consider the linear system of curves TCd-2 of degree (d — 2) 

given by the divisor J2i=i r %Qi- That is, Qi is a base point of (exact) multiplicity r j — 1 

of the linear system. Afterwards, one computes (d — 3) e-simple affine points on C h 

(see below for details), and determines the linear subsystem 7~C d __ 2 °f ^d-2 obtained 

by intersecting 7i d -2 with the linear system of (d — 2)-degree curves generated by the 

divisor $^i=i P%\ sa Y ^hat Pj := (p^i : '■ 1). If Pi, Qj would be exact points and 

singularities, respectively, of C h , then dim(7ij_ 2 ) = 1 (see Chap. 4 in [23]). However, 

in our case, since we are working with e-points we can only ensure that dim(7i d _ 2 ) > 1 

(see Theorem 2.56 in [M])- If this dimension is strictly bigger than 1, we can either take 

more e-simple points till dimension 1 is reached, or we can take an small perturbation 

of the e-points such that the effective divisor J2t=i r iQi + ^f=i Pi is m general position 

(see page 49 in [23]), and hence the dimension is 1. So, we can assume w.l.o.g. that 

dim(H d _ 2 ) = 1- Let, then, H d _ 2 (t, x, y, z) be the defining homogeneous polynomial of 
— * 

At this point, if Pj, Qj would be exact points and singularities, respectively, of 
C, the symbolic algorithm presented in Section [2] would output the parametrization 
^) = (?f,?§), where 

qi{t)x-p!{t) 



Res y (H d _ 2 (t, x, y, 1), f(x, y)) 
Res x (H* d _ 2 (t, x, y, 1), f(x, y)) 



However, in our case, Pi, Qj are not exact points, but e-points. So these rational 
functions are not, in general, polynomials. Nevertheless, considering if necessary a 
small perturbation of H d _ 2 , the quotient of the division of each numerator by its 
denominator is linear as polynomial in either x or y. Then, the idea is to determine 
the parametrization from these linear quotients. For this purpose, we will consider 
(if necessary) two perturbations, both affecting H d _ 2 . The first one will ensure that 
the degree in the resultants is the expected one, namely d{d — 2). The second will 
guarantee that the output is indeed a parametrization; i.e. that not both components 
are constants. Note that, in the exact case, these two facts are provided by the theory. 

More precisely, let H* d _ 2 {t,x,y,z) = H\(x,y,z) + t H 2 [x,y, z), and let T>i be the 
projective curve defined by Hi, i = 1,2. We recall that (1 : : 0), (0 : 1 : 0) ^ C h . 
Now, we need to ensure that either C h ,T>\ or C h ,T> 2 do not have common points at 
infinity. If this is not the case, let {Pi, . . . ,P m } be the points of C at infinity and 
K(pi, P2, x, y, z) = pix d ~ 2 + P2y d ~ 2 , where pi are parameters. Then, we consider in C 2 
the union C of the affine lines defined by P 2 (P«) + K(.Pi, P2, Pi) — 0, for i — 1, . . . , m. 
Note that, since Pj are points at infinity, the polynomials P 2 (Pj) + P(pi, P2, Ri) £ 
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C[pi, p 2 ] are not constant, and hence define lines. So, taking values for p\,p 2 (say, 
small real numbers) we consider an small perturbation that ensures that the above 
requirement is satisfied. 

Thus, in what follows we assume that T> 2 and C h do not have common points at 
infinity. Therefore, if F is the homogenization of /, by Lemma 3.1 in [I], one has that 

deg x (Res y (H d _ 2 , F)) = deg y (Res x (H* d _ 2 , F)) = d(d - 2). 

Moreover, since 7~C d _ 2 an d C h do not have common points at infinity, it holds that 

deg x (Res y (H d _ 2 (t, x, y, 1), /)) = deg y (Res x (H* d _ 2 (t, x, y, 1), /)) = d(d - 2). 

Now, we consider the polynomials 

s d—3 s d—3 

A 1 (x) = l[(x - q^-V l[(x - Phl ), A 2 (y) = - ]J(y - Pi, 2 ). 

i=l i=l i=l i=l 

Since C is e-rational, it holds that 

deg^Axix)) = deg y (A 2 (y)) = d(d - 2) - 1. 

Let Bi(x, t) := qi{t)x— Pi(t) be the quotient of S±(x, t) := Res y (H* d _ 2 (t, x, y, 1), f(x, y)) 
and Ai(x). Similarly let B 2 (y,t) := q 2 {t)x — p 2 (t) be the quotient of S 2 (y,t) : = 
Res x (H d _ 2 (t, x, y, 1), f(x, y)) and A 2 {y). Then, we output 



V(t) 



Pi(t) Pz(t) 
Qi(tYq 2 {t) 



as approximate parametrization of C. 

Intuitively one sees that, in practice, V(t) will be always a parametrization. In 
order to prove this claim, we repeat the reasoning but introducing a new perturbation 
of T~C d _ 2 . More precisely, let A = (Si, . . . , Sq) be a family of perturbing parameters and 
let 

G(A, x, y, z) = 5 iy d ~ 2 + 5 2 y d ~ 3 z + 5 3 x d ~ 2 + b±x d ~ z z + 5 5 x d ~ 3 y + 5 6 xy d ~ 3 . 

If d = 3 we take A = (Si, S 2 , S3) and G = S\y + S 2 z + <5 3 x. Observe also that in [TT] . 
Lemma 1, it is proved that for e-monomial curves, and hence for d — 3, V(t) is always 
a parametrization. Then we consider H (A, t, x, y, z) = H d _ 2 (t, x, y, z) + G(A, x, y, z); 
that is 

H**(A,t,x,y,z) = H ± (x,y,z) +tH 2 (x,y,z) + G(A,x,y,z). 

Note that we are perturbing Hi and hence H 2 keeps the required conditions on the 
point at infinity of C. 
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In this situation, repeating the above process with H and F, instead of with H d _ 2 
and F, we introduce Sf, S 2 A , Bf , R*, and P A (A, t). So 

Sf{A,x,t) = Res y (H**(A,t,x,y,l),f), S A (A, y, t) = Bes x (H**(A, t, x, y, 1), /), 

and i? A (A, x, t), i? A (A, x, t) are the quotient and the remainder of the division of S* A by 
Ax(x), respectively. Similarly, for B 2 (A, y, t), -R A (A, y, t) using S 2 and A 2 (y). Finally, 

the components of V (A,t) are the roots of 5 A (A,x, t) and B 2 (A,y,t) as univariate 
polynomials over C[A,t]. 

We start with some lemmas. 

Lemma 4.1. The leading coefficient w.r.t. x o/i? A (A,x, t) and the leading coefficient 
of B^(A,y,t) w.r.t. y, as polynomials in C(A)[£], are the same up to multiplication by 
non-zero constants in C. Furthermore, the roots are 

f Hx(a,b, 0) + G{A,a,b, 0) 
{ H 2 (a,b,0) 

Proof. Let 5 A (A, x, t) = 9l (A, i)z - Pl (A, *), and 5 A (A, y, t) = q 2 (A, t)y - p 2 (A, t). 
By hypothesis F(l, 0, 0) ^ 0, F(0, 1, 0) ^ 0. So, the leading coefficient of F w.r.t. y is 
a non-zero constant; similarly w.r.t. x. Thus, by well known properties on resultants 
(see, e.g. Lemma 4.3.1. in [23]), it holds that up to multiplication by a non-zero 
element in C: 

Res y (H**(A, t, x, y, 0), F(x, y, 0)) = (S A ) H (A, x, 0, t), 

Res x (iT*(A, t, x, y, 0), F(x, y, 0)) = (S A )"( A, y, 0, t), 

where (S^) H denotes the homogenization of S" A as polynomials in C(A, t)[x, y}. Now, 
observe that 

(S?) H (A, x, 0, t) = Ql (A, t)x^" 2 ) , (5 A ) H ( A, y, 0, t) = g 2 (A, t)y d ^ . 
Moreover, let F(x, y, 0) factor as 

d 

F(x,y,0) = Y[((3iX - any). 

i=l 

Since F(0, 1,0) ^ then aij 7^ for all i. Hence, up to multiplication by non-zero 
constants 

d 

Res y (H**(A,t,x,y, 0), F(x, y, 0)) = JJ Res y (H**(A, t, x, y, 0), fax - any) = 




13 



= (_i)^2) x d(d-2) f[iT{A, t, a,, fa, 0). 

i=l 

Analogously, 

d 

Res x (H**(A, t, x, y, 0), F(x, y, 0)) = (-1)^)^3) "Q ^" (A; ^ ^ ^ 0) _ 
So, up to multiplication by non-zero constants 

qi (A,t) = q 2 {A,t) = l\liH**(A,t, ai ,fa,0) = n!i(#iKA,0) + G(A,a 4 ,A,0) + 
tH 2 ( ai ,fa,0)). □ 

Lemma 4.2. For all A G C 6 ; deg t (5f (A , x, £)) = d and deg t (B^(A ,y,t)) = d. 

Proof. First note that deg t (B^) < d and deg t (.B A ) < d. The equality follows from 
the last equality in the proof of Lemma I4~T1 and using that H 2 (ai, fa, 0) 7^ for all i. □ 

Lemma 4.3. There exists a non-empty Zariski open subset Q ofC G such that if A G Q 
then B^(A ,x,t) and B^(A ,y,t) are primitive w.r.t. x andy, respectively. 

Proof. We assume that d > 3; if d = 3 the reasoning is analogous. Let us assume that 

Bf(A,x,t) = D(A,t)L{A,x,t), 

with deg f (.D) > 0. Then, 

Sf (A, x, t) = D(A, t)L(A, x, t)Ai(x) + (A, x, t). 

By Lemma 14.11 we know how the roots of D(A,t) G C(A)[t] are. Now for each root 
t of D (say that t is defined by P := (a : b : 0) G C ft ), deg x (5 1 A (A, x, t )) = 
deg x (i?f (A, x, to)) < d(d — 2) — 2. Let P(£o) be the projective curve defined by 
H (A, t , x, y, z) over the algebraic closure F of C(A). Then, V(t ) and C h intersect 
at infinity at an additional point different from P, or the multiplicity of intersection of 
both curves at P is at least two. We analyze each case. But first we introduce some 
additional notation. We express F, Hi, and H 2 as 

F(x, y, z) = f (x, y) + fi(x, y)z H h f d (x, y)z d , 

H(x, y, z) = h i>0 (x, y) + h it i(x, y)z^ h h i4 _ 2 (x, y)z d ~ 2 , 

where fj, hij are homogeneous of degree d—j and (d—2)—j, respectively. Moreover, we 
denote by F x , Hf, H**' x , fj, hfj the corresponding partial derivative w.r.t. x; similarly 
w.r.t. y and z. 

Let us assume that Q G T>(t ), with Q = (n : m : 0) 7^ P. This is equivalent to 

61C1 + 5 3 C 3 + 5 5 C 5 + 5 6 C 6 = C , 
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where 

Ci = b d ~ 2 h 2fi (Q) - m d - 2 h 2 , (P), C 3 = a d ~ 2 h 2fi (Q) - n d ~ 2 h 2fi (P), 

C 5 = a d - 3 bh 2t0 {Q) - n d ~ 3 mh 2t0 (P), C 6 = ab d ~ 3 h 2fi {Q) - nm d " 3 h 2fi (P), 
C = h lfi (P)h 2fi (Q) - h 2fl (P)h lfi (Q). 

Observe that h 2: o(Q) ^ 0, h 2 ,o(P) 7^ 0. Let us see that all Ci,i > 0, can not vanish 
simultaneously. Let C\ = C% = C5 = Cq = 0. We assume that a ^ 0. If a = then 
b ^ 0, and the reasoning is similar. From C3 = one has that n^0. So 

P = ( a : b : 0) = (a d ~ 2 : &a d ~ 3 : 0) = (a d ~ 2 h 2fi (Q) : 6a d - 3 /i 2 , (Q) : 0) = 

= (n d ~ 2 h 2 , (P) : n d - 3 mh 2t0 (P) : 0) = (n : m : 0) = Q, 

which is a contradiction. Therefore, if V± is the hyperplane in C 6 defined by 8\C\ + 
+ ^sCs + 5qCq = Co, for all A in C 6 \ V\ this case does not happen. 
Let us assume that the multiplicity of intersection of T>(t ) and C h at P is at least 
two. Since C h does not have singularities at infinity, this implies that both curves have 
the same tangent at P. This is equivalent to demand 

(F*(P) : F»(P) : F*(P)) = (H**' X (A,t ,P) : H**' V (A,t , P) : iT^A, t , P)). 

By hypothesis a& 7^ 0. So, by Euler's formula and taking into account that P is at 
infinity, the condition is equivalent to 

F*(P)H** >z (A,t ,P) = F z (P)H**> x (A,t ,P). 

That is equivalent to 

8 l C 1 + 5 2 C 2 + 5 3 C 3 + 5 4 C 4 + 5 5 C 5 + 5 6 C 6 = C , 

where 

a = 6 rf - 2 (A(p)^ (p) - fs(p)h 2A (p)), 

c 2 = b d - 3 fs(P)h 2fi (P), 

C 3 = a d - 3 (a/ 1 (P)/if i0 (P) - afS(P)h 2>1 (P) - (d - 2)/ 1 (P)/i 2 , (P)), 
C 4 = a d - 3 fS(P)h 2i0 (P), 

C 5 = a d -%-fS{P)h xl {P)a + / 1 (P)/if )0 (P)a - (d - 3)/ 1 (P)/i 2 , (P)), 
C 6 = b d -*{h{P)hl {P)a - h(P)h 2>0 (P) - f*(P)h 2A (P)a), 

C = fS(P)(h 2 AP)hi, (P) - h hl (P)h 2fi (P)) + h{P){hl >Q {P)h 2fi {P) - hl (P)h lfi (P)). 

Let us see that Ci,i > 0, cannot vanish simultaneously. Let G\ — ■ ■ ■ = Cq — 0. Since 
^2,0 (P) 7^ 0, and a6 7^ 0, one has that 

Ci = fx(P)hl (P) = fZ(P)h 2tl (P), 

C 2 = 0^ / C (P) = 0, 

C 6 = 0^ / 1 (P)(^ i0 (P)a-/i 2i0 (P))=0. 
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Note that /i(-P) 7^ 0, since otherwise it would imply that F Z (P) = and using that 
F X (P) = f*(P) = and that (by Euler's formula) F y (P) = 0, one would deduce that P 
is a singularity of C h which is excluded by hypothesis. So, the first and second equalities 
imply that h^^P) = and this yields to (using the last equality) h 2 p(P) = which 
is a contradiction. Therefore, out of the hyperplane defined in C 6 by X^=i = Co, 
this case cannot happen. 

For each point of C at infinity we generate the hyperplanes described above and 
corresponding to each one of the two cases. Let V be the union of all of them, and 
let fli = C 6 \ V. Repeating the same reasoning with (note that G is symmetric in 
terms of x and y), we get Q 2 . Finally, let O = Qi n 02- □ 

Now, the next theorem follows directly. 

Theorem 4.4. There exists a non-empty Zariski open subset Q of C 6 such that if 
A Gfi then V (A ,t) is a rational parametrization of a rational curve of degree at 
most d. 

Proof. Taking Q as in Lemma l4~3l we ensure that P A (Ao, t) is a rational parametriza- 
tion. By Lemmas 14.11 and 14.21 we get that the degree of the curve is at most d. □ 

Remark 4.5. Let 7iX be the linear system of (d — 2)-degree curves defined by 
H**(A , t, x, y, z). If no perturbation is needed, i.e. A = 0, then Ti.$* = 7~C,- 2 , and 
hence it is generated by the effective (exact) divisor Yli=i r^Qi + Ylt=i Pi- Now, if we 
identify (as usual) C 2 with M 4 and we consider the perturbing parameters Si as real 
variables, it holds that for each Qi (similarly of Pi) there exists p(Qi) > such that for 
almost every element (Q*, A ) in the open Euclidean disk of M 10 , of center (Qi, 0) and 
radius p(Qi), V(A ,t) is a parametrization and Q* is an e-point of e-multiplicity (at 
least) Ti of the curves defined by H 1 (x, y, z) + G(A , x, y, z) and H 2 (x, y, z) (recall that 
H = H 1 + G + tH 2 ); i.e. of the generating curves of H A . This can be seen by applying 
Theorem 14.41 taking 5i small enough to ensure that = \\Hi + G(A , x, y,z)\\, and 
noting that if M is any of the derivatives of Hi +G and H 2 involved in the e-multiplicity, 
then 1 1 JW 1 1 § is a continuous function that vanishes at (Qi,0). □ 

Finally, and before outlining the algorithm, we briefly describe how to proceed with 
the selection and computation of the (afline simple) e-points Pi. We first observe that, 
in general, an e-point can be computed by solving {f(x,y) = 0, ax + [3y = p}, where 
a,j3,p G C, under fixed precision e||/||. However, we are intersected in working with 
either real e-points or pairs of conjugate complex points. We can always compute all 
points, but at most one, in pairs of conjugate complex points. For choosing real points 
one can always analyze the roots of the discriminant of / (see Theorem 7.7 in |24] ) . 
On the other hand we have observed, in our examples, that taking (when possible) 
the simple e-points as (afline) e-ramification points (see Def. 11.31) the error distance 
between the original curve and the output curve decreases. So we tend to use first such 
points. Finally, one has to take care of the fact that a chosen e-point can be too close 
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(i.e. in the same cluster) to an e-singularity or to a previously computed e-point, and 
hence identifiable with it. To avoid this, whenever a new simple e-point is computed 
we check whether it belongs to the cluster of the others points. 

The above process provides the following approximate parametrization algorithm 
for deciding whether a real e-irreducible (with proper degree) plane algebraic curve C 
is e-rational, and in the affirmative case, compute an approximate parametrization. 
Recall that we assume that C does not have exact singularities at infinity, and that 
(0 : 1 : 0),(1 : : 0) £ C h . If this last condition fails, one may consider an afline 
orthogonal change of coordinates to achieve the requirement. 

Approximate Parametrization Algorithm 

• Given a tolerance e > and an e-irreducible polynomial f(x,y) G R[x, y], of 
proper degree d > 2 (for d — 1 it is trivial, if d = 2 apply [17]), without exact 
singularities at infinity, not passing through (0 : 1 : 0), (1 : : 0), and defining a 
real plane algebraic curve C; let F(x,y, z) be the homogenization of /. 

• Decide whether C is e-rational and in the affirmative case 

• Compute a rational parametrization V(t) of a curve C close to C. 

(1) Compute the cluster decomposition {<C/Mster n (Qi)}j=i,..., s of Sing e (C); say Qi = 
(<?i,i : <?i,2 : 1)- 

(2) If £i=i r i(Vi - 1) ^ (d - l)(d - 2), RETURN "C is not (affine) e- rational 1 . If 
s = 1 one may apply the algorithm in |17j . 

(3) Determine the linear system Hd-2 of degree (d — 2) given by the divisor Yli=i r iQi- 

(4) Compute (d — 3) e-ramification points {-Pj}i<j<<j-3 of C; if there are not enough 
e-ramification points, complete with simple e-point. Take the points over R, or 
as conjugate complex points. After each point computation check that it is not 
in the cluster of the others (including the clusters of Qi); if this fails take a new 
one. Say P t = (p t>1 : p ij2 : 1). 

(5) Determine the linear subsystem 7i* d _ 2 of Tid-2 given by the divisor Y^t=i -^*- Let 
H*(t, x, y, z) = Hi(x, y, z) + tH 2 (x, y, z) be its defining polynomial. 

(6) lf[gcd(F(x,y,0),H 1 (x,y,0)) ^ 1] and [gcd(F(x, y, 0), H 2 {x, y, 0)) + 1] replace^ 
by H 2 + pix d ~ 2 + p 2 y d ~ 2 , where p±, p 2 are real and strictly smaller than e. Say 
that gcd(F(x, y, 0), H 2 (x, y, 0)) = 1; similarly in the other case. 

(7) Set <&i = • • • = 8 6 = 0. 

(8) If d > 3 then g := 5iy d ~ 2 + 5 2 y d ~ z + 5 3 x d ~ 2 + 5 4 x d - 3 + 5 5 x d -' 3 y + 5 6 xy d - 3 else 
£ := 8 x y + 5 2 + 5 3 x. 
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(9) Sfat) = Res y (H*(x,y,l) + g,f) and S 2 (y,t) = Res x (H*(x,y, l)+g,f). 

(10) a ± (x) = nu(^-ft,i) r4(Pi - 1) n£f(^-ft,i) > 
^2(1/) = nuG/ - ^r (r '- 1] xitliy-Pi?)- 

(11) For 2 = 1,2 compute the quotient 5j of Si by w.r.t. either x or y. 

(12) If the content of B\ w.r.t x or the content of B2 w.r.t. y does depend on t, take 

. . . , 8$} as small real numbers (strictly smaller than e) and go to Step 8. 

(13) Determine the root Pi(t) of B\ as a polynomial in x and the root p 2 (t) of B2 as a 
polynomial in y. 

(14) RETURN V(t) = (pi(t),p 2 (t)). 

The next theorem states the main properties of the curve output by the algorithm. 
But first, we need the following technical lemma. 

Lemma 4.6. Let L be the algebraic closure of C(t), and €1,62 two plane projective 
curves over L with defining polynomials Gi(x,y, z),G2(x,y, z) G C[t][x, y, z], respec- 
tively. If there exist K,W,Le C[t] [x, y, z] such that KG\ + WG 2 = zL, and 

(1) G 1 (x J y J 0)G 2 (x,y,0)^0, 

(2) gcd(G 1 (x,y,0),G 2 (x,y,0)) = l, 

then either z divides K and W or there exist U\, U 2 , U3 G C [£][#, y, z] such that 

L = U x G x (x, y, 0) + U 2 G 2 {x, y, 0) + zU 3 . 

Proof. If z divides K, then z divides WG 2 , and by (2) z divides W. So let us assume 
that z does not divides K, and let us denote by G° the polynomial Gi(x, y, 0); similarly 
with K°, W°. Then, K Q G\ + W G° 2 = 0. Since G° ^ and gcd(G°, G° 2 ) = 1, then G\ 
divides W° and G° 2 divides K°. Let K° = AiG°, W° = A 2 G° V So (Ai + A 2 )G?G 2 I = 0, 
and since G® 7^ 0, one gets Ai + A 2 = 0. Now, we write 

K = K° + zK, W = W° + zW, G i = G° i + zG u 

where K,W,G t G C[t][x,y,z\. Then, KW l + WG 2 = z^Ji +G° 2 W +z(KG l + WG2))■ 
u 

Theorem 4.7. The rational curve C, output by the algorithm, and C have the same 
points at infinity, and deg(C) < deg(C). 
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Proof. The fact on the degree follows from Theorem 14.41 For the rea- 
soning of the rest of the proof, we can assume w.l.o.g. that no perturba- 
tion A is required (i.e. A = 0) in the execution of the algorithm. Let 
H* d _ 2 (t, x, y, z), F(x, y, z), S x (x, t), S 2 (y, t), A x {x), A 2 (y), B x (x, t), B 2 (y, t), Ri(x, t), and 
R 2 (y,t) be defined as above. Let B x := q x (t)x —pi(t), B 2 := q 2 (t)y — p 2 (t), and recall 
that Ri is the remainder of the division of Si by Ai. Furthermore, by Lemma 14.11 
qi (t) = Xq 2 (t), with A G C*. By Lemma S3 deg t (5i) = deg t (B 2 ) = deg(F) = d, and, 
by Lemma SiS, gcd(q 1: p x ) = gcd(q 2 ,p 2 ) = I. So, 

V H (t) := (A-V(t) : P2(t) : fe(t)) 

parametrizes the projective closure of C. Furthermore, since deg(pj) < deg(g 2 ), then 
all points of C at infinity are reachable by V (t) (see [23]). In addition, we note that 

deg {x>y>zy (H* d _ 2 ) = d - 2, deg(^) = d(d - 2) - 1, deg fey} (i2,-) < d(d - 2) - 2. 

Moreover, if m H (x,y,z,w) denotes the homogenization of m(x,y,w) as a polynomial 
in C[u>][a;, y], we have that 

S*(x,z,t) = Res y (H* d _ 2 (t,x } y,z) : F(x } y,z)) = B H (x, z, t)A? (x, z) + R? (x, z, t)z n \ 

S 2 {y, z, t) = Res x (H d _ 2 {t, x, y, z), F(x, y, z)) = B*(y, z, t)A*(y, z) + R*(y, z, t)z n \ 

where rij + deg(Rf) = d(d — 2), j = 1, 2. So rij > 2. Also, we denote by and 
the set of points at infinity C and C respectively. By resultant properties, there exist 
polynomials M x , N x , M 2 , N 2 G C[t, x, y, z] such that 

Mjt^ + NiF^sf, z=l,2. 

So, 

yA*S? - XxAfS 2 = zA?A%(\xp 2 - y Pl ) + z n *R 3 , 

where n 3 > 2 and R3 a polynomial; namely z n3 R 3 = yA 2 z ni Ri — \xA x z n2 R 2 . On the 
other hand, if K = yA%M 1 - XxA\ H M 2 and, W = yA 2 H N 1 - \xA?N 2 , then 

yA%S? - XxAfSi? = K(x, y, z, t)(H x + tH 2 ) + W(x, y, z, t)F. 

Therefore, z divides the right hand side of the above equation. We now check that H x + 
tH 2 and F satisfy the hypothesis of Lemma 14.61 Since F is irreducible and non- linear, 
F(x, y, 0) ^ 0. Moreover, if H x (x, y, 0) + tH 2 (x, y,0) = then H 2 (x, y, 0) = and this 
implies that T> 2 contains all the points at infinity of C h , which is a contradiction. Finally, 
if gcd^z, y, 0)+tH 2 (x, y, 0), F(x, y, 0)) ^ 1, then gcd(# 2 (x, y, 0), F(x, y, 0)) ^ 1, and 
this implies that T> 2 and C h share points at infinity. Therefore, applying Lemma 14.61 
one deduces that either there exist polynomials M 3 ,N 3 G C[£][x, y, z] such that 

M 3 H d _ 2 + N 3 F = A^Af(Xxp 2 - yp x ) + z n *R 3 , 
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where > 1, or there exist polynomials U±, U 2 , U3 G C[£][x, y, z] such that 

Uj? d _ 2 (t, x, y, 0) + U 2 F(x, y, 0) + zU 3 = A?A 2 H (\xp 2 - yp 1 ) + z n *R 3 . 

In this situation, using C V (C), we first observe that Card(Coo) is less or equal 
to the number of different roots of q 2 (t) and, by Lemma 14.11 this number is less or 
equal to Card(Coo). So, Card(Coo) < Card(Coo). Now, we prove that C Coo, from 
where one concludes the proof. Let P = (x : y : 0) G C^, and let t be the root of q 2 
generated by P (see Lemma H~Tj) . So, H d _ 2 (t ,x ,y ,0) = F(x ,yo,0) = 0. Applying 
the corresponding equality above, and using that n 4 > 1, we get 

A? (x , 0)A% (y , 0)(Xx p 2 (t ) - y oPl (t )) = 0. 

Moreover, since (1 : : 0), (0 : 1 : 0) ^ C h then xoyo ^ 0, and hence 
Ay (x , 0)A 2 (y , 0) ^ 0. So, Xx Q p 2 (t ) = yoPi(t ). In addition, pi(t )p 2 (t ) ^ be- 
cause gcd(g 2 ,pi) = 1 = gcd(g2,P2)- Therefore, 

V H (t ) = (\- 1 Pl (t ):p 2 (t ):0) = 

= (z/oA~V(*o) : yoP2(t ) : 0) = (x p 2 (t ) : y p 2 (t ) : 0) = (z : y : 0) = P. 

□ 



5 Displaying Examples. 



In this section we present several examples (the degrees are 5,6,7) to illustrate the 
algorithm. These examples have been computed in Maple. 

0.01 and C the curve of proper degree 5 defined by the poly- 



Example 5.1. Let e 

nomial (see Figj3]): 

f(x,y) = 0.006521014507x 4 + 0.006521014507a; V - 0.3174429862a; 3 
0.006521014507a; 4 y + 0.03536521618?/ 4 + 0.008903520149xV - 0.1541837293y 3 
0.3561209555a; 2 - 0.2351465855?/ 2 + 0.01517989182a;y 4 + 0.006177658243?/ 5 
0.006521014507a; 5 - 0.6503293396xy + 0.006521014507a; 3 y - 0.6965951291xy 2 
0.1751383118a;V + 0.1487535027a;?/ 3 - l.x 2 y + 0.0000658688334. 

First we compute the e-singularities of C, obtaining the e-singular locus Sing e (C) 
Si U S 2 U S 3 : 



51 = {Qi 

Q3 
Qa 

Q 5 
Q7 

5 2 = {Q S 

5 3 = 0. 



-3.999854219,2.000094837), Q 2 = (0.,0.), 
0.9998153818, -2.999388343), 

-2.001190360 + 0.05414244305?:, 3.001898191 - 0.08039416354i), 
-1.980207988,3.002780607), Q 6 = (-2.019931003,2.997118979), 
-2.001190360 - 0.05414244305z, 3.001898191 + 0.08039416354z)}, 
-2.000000001, 3.000000001)}, 



+ 

+ 
+ 



20 



o 



o 



o 



o 



o 



o 



o 







Figure 2: Left: Cluster Decomposition of the Singular Locus. Right: Cluster Decomposition 
of the Singular Locus with two e-ramification points. 

Moreover, the cluster decomposition of the singular locus is (see Fig. Left): 

<£luster 2 {Qi) = {Qi}, <£luster2(Q2) = {Q2}, £luster 2 {Q3) = {Q3} and 

Cluster 2,{Qs) = {Qa, Q5, Q7, Qs}- 
We observe that C is e-rational. Following Step 4 in the algorithm we ob- 
tain two e-ramification points, namely P\ = (3.437938023,4.260660564), P 2 = 
(7.712891931, 1.573609575). We note that these points are not in the cluster of each 
other and they are not in the clusters of the cluster decomposition of the singular locus 
(see Fig. El Right). 

Finally, the algorithm outputs the parametrization V(t) = (^^, "where (see Fig. 
[3] to compare the input and the output curves): 

p x (t) = -0.5918689071 10" 29 (0.7256428750 10 579 t + 0.1009796140 10 581 i 4 + 
0.4757134093 10 580 t 5 + 0.3531628351 10 580 £ 2 + 0.8491037424 10 580 t 3 + 
0.5883163866 10 578 ) 

p 2 (t) = 0.3851669500 10~ 31 (0.1621127956 10 583 t 3 + 0.1491645111 10 582 t + 
0.6997743856 10 582 t 2 + 0.8444468165 10 582 t 5 + 0.1252710479 10 581 + 
0.1858263849 10 583 t 4 ) 

q(t) = 0.1265532998 10 551 t 3 + 0.13722 1 7100 10 551 t 4 + 0.1260572385 10 549 + 
0.1356321818 10 550 t + 0.5851539780 10 550 * 2 + 0.5967959572 10 550 i 5 

We note that the algorithm did no require perturbing H d _ 2 - □ 

Example 5.2. Let e = 0.004 and C the curve of proper degree 6 defined by the 
polynomial (see Figj5]): 
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Figure 3: Left: Input (in dots) and output curve in Example 15.11 Right: A zoom at 
the origin 



f(x, y ) = -0.5499999998 x + 0.9999999999 y + 0.00002677376171 + 0.006666666664 x 2 y 3 + 
0.006666666664 x 4 y + 0.3799999999 x y 2 - 0.4133333332 x y 3 + 0.006666666664 x 5 y - 



0.01999999999 xy b 



0. 1066666667 x A y - 0.07000000000 x^ y 4 + 0.8066934397 x y 



0.03333333332 x 2 y + 0.03999999998 x 3 y 3 + 0.5466666665 x 2 y 2 + 0.1133333333 x 3 y 2 
0.04999999998 xy 4 - 0.5333333332 y 3 + 0.006666666664 x 6 + 0.006666666664 x A y 2 
0.006666666664 y 6 - 0.6700000000 x 2 - 0.1766666666 x 4 + 0.3599999999 y 4 
0.4699999998 x 3 - 0.006666666664 x 5 - 0.6066666665 y 2 - 0.03999999998 y 5 . 

We get the e-singular locus Sing e (C) = Si U S2 U S3 where 



Si 



s 2 

S3 



{Qi- 
Q2-- 
Q3-- 

Qa~- 
Qe~- 
Q7-- 

Q 8 ~- 

Q9-- 

Qw 
Qn 
Qn 

{Ql4 



-1.994232333, 
-2.000005299 ■ 
-2.000014217- 
-2.003547061, 
-1.996418580, 
-2.000014217- 
-2.000005299 ■ 
1.000036272 + 
(5.999999669, - 
(1.000036272 - 
(0.9978910941, 
(-2.000000001 



1.005043048), 

f 0.005645280797i, -1.000026945 - 0.0002822677587i), 
f 0.004619269427i, 1.000004775 - 0.003559494332i), 
-1.006293429), Q 5 = (-2.005740475,0.9948974977), 
-0.9936748962), 

- .004619269427^, 1.000004775 + 0.003559494332i), 

- 0.005645280797z, -1.000026945 + 0.002822677587*), 
0.008596901071i, 2.000017052 - 0.003059926359z), 
-2.999998564), 

0.008596901071z, 2.000017052 - 0.003059926359z), 
1.994329680), Q 13 = (1.002094534,2.005650021)}, 
, 1.), Q 1S = (-2., -1.000000005), Q w = (1, 2.)}, 



The singular cluster decomposition is (see Fig. HI Left): 

£luster 2 (Qw) = {Qw}, <tluster 3 (Qu) = {Qi, Q 3 , Q 5 , Q 7 , Q u }, 

<tluster 3 (Qi 5 ) = {Q 2 , Qi, Qe, Qs, Q15}, and £luster 3 (QiQ) = {Q 9 , Qn,Qi 2 , Q13, Qw}- 
We observe that C is e-rational. In Step 4 we obtain three e-ramification points: 
Pi = (-1.330235522, 0.9268173641), P 2 = (-1.979908167, 0.02661222172), and P 3 = 



22 



(-2.700785807, -0.07757312293). We note that these points are not in the cluster of 
each other and they are not in the clusters of the e-singularities (see Fig. [4j Right). 
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Figure 4: Left: Cluster Decomposition of the Singular Locus. Right: Cluster Decomposition 
of the Singular Locus with two e-ramification points. 



The algorithm outputs the parametrization V(t) = (^r, *W£)) where (see Fig. [5] to 
compare the input and the output curves): 
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Figure 5: Left: Input (in dots) and output curve in Example I5.2I Right: A zoom at 
(-2,-1) 



p x (t) = -0.2992985374 10~ 13 (-0.4498780650 10 665 t 2 - 0.1104625259 10 666 t + 

0.3823432112 10 663 1 5 - 0.4349945664 10 664 1 3 - 0.89 77530140 10 665 + 
0.2445532363 10 662 t 6 + 0.1487140379 10 664 t A ) 

p 2 (t) = -0.5410017657 10" 14 (0.4254697372 10 662 t 6 - 0.4450231957 10 665 t 3 - 

0.8325944623 10 665 1 2 + 0.3137371087 10 665 1 - 0.1036628422 10 663 1 5 - 
0.7148994612 10 664 t A + 0.1592853294 10 666 ) 



23 



q(t) = 0.3388387927 10 649 t 5 
0.4492917291 10 651 1 + 
0.4139801407 10 651 . 



+ 0.1633128101 10 651 t 2 + 0.3082510569 10 648 t 6 + 
+ 0.3531547270 10 650 1 3 + 



0.1270749205 10 650 t A 



We note that the algorithm did no require perturbing H d _ 2 . 



□ 



0.001 and the curve C of proper degree 7 defined 



Example 5.3. Let us consider e 
by the polynomial (see Figj6]): 

f(x,y) = 0.005242164122x + 0.0000006109092905 + 0.4234041949y - 
0.1626221914xV - 0.006150324474x 5 y 2 - 0.009115378696y 3 + 
0.1726592957y 4 - 0.005178781717x 5 + 0.0006102983812y 6 + 



0.007271579029x 4 
0.03335531981xy 3 
0.06359239287x 2 y 
0. 1735938027a; V + 



0.005178781717x° + 
- 0.009049345878x 3 y 4 
+ 0.07030423460xV 
3 + 0.0006102983812xy 5 
0.009673386920xy 



0.05219720755xy 4 
0.01468749412a; 3 
0.7056394692y 2 



0.01517536020xV 
0.02396026447x 5 y 
0.05961614786x 4 y 2 
0.3997312415x 2 i/ 

0.01504641433a; V - 0.0002034043985x Y + 0.0007152781730x tl y + 0.009647777670xy 6 
0.01996027092x 2 - 0.001858780227x 6 - 0.008636725103x 2 y 5 - 0.002554427076y 7 . 

The e-singular locus is Sing e (C) = Si U S2 U S3 where 



+ 0.02810421594y 5 + 
+ 0.9999999999xy 2 + 
+ 0.06037915453x 4 y - 
- 0.1183998159x 3 y - 

6, 



Si 



So 



{Qi- 
Q2-- 
Q3-- 

Qa~- 

q 5 = 

Q 6 = 

Q7-- 

Qs~- 

Q9-- 

Qio 

Qn 

Q12 

Ql3 
Ql4 

Q15 
Q16 
Qn 

Q18 

Q19 

{Q20 
Q21 
Q22 



0.000335729506H), 

0.000335729506H), 
0.0003352187676i), 



4.998181206 + 0.00046392220802, 6.997094116 
4.998181206 - 0.0004639222080i, 6.997094116 
5.001816967 + 0.0004635470406i, 7.002902744 
0.9999998537, -3.000000118), 

5.001816967 - 0.0004635470406*, 7.002902744 + 0.0003352187676*), 
-2.000211362 + 0.00008683312445*, -0.001218550314- 0.9986341029*'), 
-2.000211362 - 0.00008683312445*, -0.001218550314 + 0.9986341029i), 
1.998594026 + 0.0001453051485*, -0.0005005279912 - 0.9994646423^), 
-1.998594026 - 0.0001453051485*, -0.0005005279912 + 0.9994646423*'), 
1.000001333, -0.00000005539458512), 

-2.001405450 - 0.0001446569643*', 0.0005000190394 - 1.000535463*), 
-2.001405450 + 0.0001446569643i, 0.0005000190394 + 1.000535463*), 
-1.999787395 - 0.00008416274464i, 0.001216340837 - 1.001365249*"), 
-1.999787395 + 0.00008416274464i, 0.001216340837 + 1.001365249i), 
4.997608917 - 0.001980994691*, 1.999734804 - 0.002346469999*'), 
4.997608917 + 0.001980994691i, 1.999734804 + 0.002346469999*), 
-3.999997183, 1.999998082), 

5.002393988 - 0.001973712849*, 2.000267815 - 0.002341789270*'), 
5.002393988 + 0.001973712849^, 2.000267815 + 0.002341789270*')}, 
-2.000000398 - 0.0000003109941563*, 
-0.0000005243819124 - 0.9999997083*), 
-2.000000398 + 0.0000003109941563i, 
-0.0000005243819124 + 0.9999997083i), 

5.000000495,2.000000179), Q 23 = (4.999999480,6.999999337)}, 
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The cluster decomposition of the singular locus is (see Fig. HI Left): 

€luster2{Qi) = {Qa}, <£luster 2 (Qio) = {Qio}, Cluster 2 (Q '17) = {Qn}, 
Clustery (Q23) = {Qi, Q2, Q3, Q5, Q23}, £luster 3 (Q 20 ) = {Q 6 , Q 8 , Qn, Qi8, Q20}, 
Cluster 3 (Q 21) = {Q7, Q9, Q12, Qu, Q21}, and £luster 3 (Q 22 ) = {Q15, Qw, Qi8, Q19, fe}- 
We observe that C is e-rational. Now, we obtain four e-ramification points: 
Pi = (-2.972405737,-7.933174980), P 2 = (23.79950366,17.84891277), P 3 = 
(-10.06218879, 1.300686562) and P 4 = (24.47385001, 17.37936091). We note that these 
points are not in the cluster of each other and they are not in the clusters of the e- 
singularities. Finally, the algorithm outputs a parametrization (for space limitation we 
do not include it here). In Fig. [6] we plot the input and the input curve. 




Figure 6: Left: Input (in dots) and output curve in Example 15.31 Right: A zoom at 
(-4,2) □ 

5.1 Empirical Analysis of Error 

The aim of this section is to analyze empirically the performance of the algorithm 
proposed. A good performance would mean to obtain an output curve C close to the 
input curve C and by close we mean that C is contained in the offset region of C at a 
small distance and viceversa (see [17]). To estimate the distance, between the curves 
C and C, we designed the next method (see Fig. [71 Left): 

1. We randomly generate a set £ of (affine) e-points on the input curve C as follows. 
Fix real numbers a, b, and take n random integer values G [a, b], i — 1, . . . , n. 
Let £\ be the set of intersection points of the curve C with the lines x = oti. We 
also take n random integer values G [a, b] and we intersect the curve C with 
the lines y = ati to obtain the set of points £2. We set £ = £\ U £2- 

2. Let r be a positive integer and r = {— | k — 1, . . . , r}. For each e-point P in £ 
and each fle8 r let Lp^ be the line through P in the direction of {cos (6), sin{9)). 
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3. For each P G £ we compute d P ^ = min{||P — Q\\ 2 \ Q G Uo e e r (Lp,e H C)}. 

4. Let T> = {dpc \ P G £}. We compute the mean value /i of the elements of T> as 
well as the statistical standard error p. We can say that the estimated distance 
is, in average, in the interval [p, — 1.96 p, p + 1.96 p\. 

Given the curve C of degree d = 4 defined by the irreducible polynomial 

f(x, y) = 1.000065?/ 2 + 1.00000028y 3 + y 4 + 1.000065sy - 11.49999972x?/ 2 

+xy 3 + 0.760065x 2 + 5.74000028a; 2 ?/ + 3.69xV - 0.75999972x 3 - 3.12x 3 y (1) 

+0.19x 4 + 0.01x + 0.01?/. 

For e = 0.01 the algorithm concluded that C is e-rational and returned a rational 
parametrization V(t) = (Pi{t) ,p 2 {t)) which corresponds to the rational curve C with 
implicit equation 

J{x,y) = 0.01642553x 4 + 0.06494377x 2 + 0.08804654?/ 2 - 0.06552116x 3 
+0.08169391y 3 + 0.091025077x?/ + 0.49976135x 2 y - 0.99999999xy 2 + 0.08645018?/ 4 
+0.319001 18xV - 0.26972458x 3 ?/ + 0.08645019xy 3 - 0.00001398 + 0.00078781x 
+0.0007408y. 

With the above notation, we take [a, b] = [—100, 100], n = 15. The number of points 
used to compute the experimental distance is \S | = 120. Set the number of lines 
going through each point equal r = 10. The mean value of the positive integers in 
T> is p = 0.007541 and the statistical standard error p = 0.000855. We can conclude 
that the estimated distance is, in average, in the interval [pi — 1.96 p,p + 1.96 p] = 
[0.005866,0.009217]. 

For the examples treated in Section [U the corresponding figures show that input 
and output curves look very close. We show below that the estimated distance between 
C and C is small compared to the value of e used in each one of the examples. For this 
purpose, we execute the previously described method to estimate the distance between 
C and C taking r = 10, and [a, b] = [—100, 100]. Let | £ \ be the number of points taken 
to compute the experimental distance. The results obtained in the examples are: 

1. Let C and C be the curves in Example 15.11 e = 0,01, | £ \= 150. Then, /i = 
0.070953, p = 0.012655. So, the estimated distance is, in average, in the interval 
[0.046149,0.095757]. 

2. Let C and C be the curves in Example 15.21 e = 0,004, | £ \= 180. Then, 
/i = 0.002425, p = 0.000100. So, the estimated distance is, in average, in the 
interval [0.002228,0.002621]. 

3. Let C and C be the curves in Example I5.3[ e = 0,001, | £ \— 210. Then, 
/i = 0.006560, p = 0.000609. So, the estimated distance is, in average, in the 
interval [0.005366,0.007754]. 
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Until now we have empirically measured the distance in the examples used in paper. 
We describe next a different experiment. We randomly generate a set of curves and for 
each curve we estimate its distance to the output curve given by our algorithm. Our 
experiments are satisfactory and allow us to think about a theoretical treatment of this 
fact as a future project. We explain next how the family of curves was constructed. 
We fix three points P 1 = (2 : : 1),P 2 = (0 : : 1) and P 3 = (1 : 1 : 1) in P 2 (c). We 
consider the linear system of curves of degree 4 defined by the divisor 2 Pi + 2P 2 + 2P3. 
Its defining polynomial is 

G(x, y, z) = u 2 y 2 z 2 + u 3 y 3 z + u 4 y 4 + u 5 xyz 2 + (-2u 2 - 3u 3 
— 4it 4 — l/2« 5 — 2u 6 )xy 2 z + u G xy 3 + uix 2 z 2 + (— 3/2u 5 + 2u 3 , . 

+4u 4 + 2u 6 - ui)x 2 yz + (u 2 + u 3 + l/2u 5 + 1/4^ + u 4 )x 2 y 2 ^ ' 

—Uix 3 z + (l/2u 5 — u 3 — 2n 4 — u 6 + l/2ui)x 3 y + l/AuiX A 

For j = 1, . . . , 6 and % = 1, . . . , 10 let r^- be a random integer number in the interval 
[0, 100]. We obtain 60 different polynomials Gy(x, y,z), j = 1, . . . , 6, i = 1, . . . , 10 
setting 



in equation ([2]). Given i £ {1, ... ,6} and j G {1, . . . , 10} we obtain a random pertur- 
bation gij(x,y) G M[x, y] of Gy(x,y, 1) as follows 



9ij(x,y) = Gi j (x,y,l) + e-^-(x + y)+e 2 -^-(x 2 + xy + y 2 ) + e 3 -^-( y x 3 + x 2 y + xy 2 + y 3 ) 



where ri, r 2 , r 3 are integer numbers taken randomly in the interval [0, 100] and e = 0.01. 
The polynomials gij(x,y), j = 1, . . . , 6, i — 1, . . . , 10 have proper degree 4 and define 
60 curves dj verifying (1 : : 0), (0 : 1 : 0) £ and without (exact) singularities at 
infinity. Using our algorithm we conclude that 28 of the 60 curves are e-rational. We 
show those curves in Fig. [TJ Right. The implicit equation of the curve Cn is equation 



With the notation used in the method described above to compute the experimental 
distance, we take [a, b] = [—100,100] and n = 15 so that the number of points used 
to compute the distance is expected to be | £ \= 120. Set the number of lines going 
through each point equal r = 10. We compute fi,p and = [// — 1.96 p, /x + 1.96 p] 
for each one of the 28 e-rational curves to obtain the next table. 
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Figure 7: Left: Illustration of the method. Right: The 28 e-rational curves Cij ran- 
domly generated . 



M 


P 




0.007541 


0.000855 


[0.005866,0.009217] 


0.006977 


0.001184 


[0.004656, 0.009299] 


0.006977 


0.001184 


[0.004656, 0.009299] 


0.003577 


0.000503 


[0.002592, 0.004563] 


0.004011 


0.000553 


[0.002928, 0.005094] 


0.006007 


0.000808 


[0.004423, 0.007590] 


0.004239 


0.000844 


[0.002584, 0.005894] 


0.005758 


0.000585 


[0.004610, 0.006905] 


0.002882 


0.000224 


[0.002442, 0.003322] 


0.005477 


0.000756 


[0.003996, 0.006958] 


0.003123 


0.000437 


[0.002266, 0.003979] 


0.004752 


0.000359 


[0.004049, 0.005455] 


0.001453 


0.000148 


[0.001163, 0.001744] 


0.004956 


0.007123 


[0.035599, 0.063522] 


0.001049 


0.000113 


[0.000827, 0.001272] 



P 


P 




0.006807 


0.000385 


[0.006051, 0.007563] 


0.100902 


0.013253 


[0.074926, 0.126879] 


0.003049 


0.000254 


[0.002551, 0.000355] 


0.003924 


0.000212 


[0.003508, 0.004339] 


0.003995 


0.000549 


[0.002919, 0.005072] 


0.008330 


0.000806 


[0.006749, 0.009911] 


0.005638 


0.000536 


[0.004586, 0.006690] 


0.003020 


0.000316 


[0.002399, 0.003639] 


0.000854 


0.000091 


[0.000677, 0.001032] 


0.004077 


0.000274 


[0.003540, 0.004614] 


0.035130 


0.005220 


[0.024898, 0.045361] 


0.006209 


0.000619 


[0.004996, 0.007423] 


0.013406 


0.001179 


[0.011094, 0.015718] 


0.009037 


0.0006872 


[0.007691, 0.010385] 
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